Universal linear Bogoliubov transformations through one-way quantum computation 
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We show explicitly how to realize an arbitrary linear unitary Bogoliubov transformation (LUBO) 
on a multi-mode quantum state through homodyne-based one-way quantum computation. Any 
LUBO can be approximated by means of a fixed, finite-sized, sufficiently squeezed Gaussian cluster 
state that allows for the implementation of beam splitters (in form of three-mode connection gates) 
and general one-mode LUBOs. In particular, we demonstrate that a linear four-mode cluster state 
is a sufficient resource for an arbitrary one-mode LUBO. Arbitrary input quantum states including 
non-Gaussian states could be efficiently attached to the cluster through quantum teleportation. 



I. INTRODUCTION 

The cluster model of quantum computation, or one- 
way quantum computation [H, 0], is an alternative ap- 
proach to the standard circuit model for quantum com- 
puting In the cluster model, a special type of en- 
tangled state is used as a resource for cluster computa- 
tion. These resource states are known as cluster states. 
A cluster computation is basically a sequence of elemen- 
tary, 'half teleportations 0, Q where quantum informa- 
tion is not only transmitted through a cluster state but 
also manipulated in any desired way depending on the 
specific choice of the measurement bases at each telepor- 
tation step. As opposed to standard-teleportation-based 
schemes, the measurements in a cluster computation are 
all local (subsequently performed on the individual nodes 
of the cluster). In order to achieve universal quantum 
computation using a fixed cluster state, active feedfor- 
ward is needed, where the measurement bases of subse- 
quent measurements have to be adjusted according to the 
outcomes of the earlier measurements. 

Cluster states and cluster computation were originally 
proposed for discrete variables (DV), namely qubits 
More recently, the cluster-state model was then extended 
to the regime of continuous variables (CV) 0,0, m 
which universal cluster states can be approximated by 
experimentally highly accessible Gaussian multi-mode 
squeezed states of sufficiently many quantized optical 
modes (qumodes). Both for DV and for CV, the cluster- 
state model is known to be equivalent to the circuit model 
in the sense that any finite-dimensional (qubits) as well 
as any infinite-dimensional (qumodes) operation can be 
efficiently realized in a cluster-based scheme. 

For DV, an arbitrary single-qubit rotation (unitary) 
can be exactly decomposed into three elementary single- 
qubit rotations [3(- Therefore, even though the whole 
set of single-qubit unitaries is continuous, concatenating 
three elementary (but continuous) single-qubit rotations 
in a three-step cluster computation using a linear four- 
qubit cluster state is sufficient to achieve universality in 
the single-qubit space. Such elementary rotations by gen- 



eral angles would include so-called non-Clifford gates; in 
this case, feedforward is required during the cluster com- 
putation. As a result, provided that the continuous, ele- 
mentary single-qubit rotations can be implemented in an 
error-resistant fashion, any multi-qubit unitary can be 
performed by connecting sufficiently many linear four- 
qubit clusters by vertical wires through which a fixed 
two-qubit entangling gate can be applied when needed. 

In the case of CV, there are various subtleties, even 
in theory. First, independent of the cluster model, an 
arbitrary single-qumode transformation (represented by 
a Hamiltonian which is an arbitrary polynomial of the 
qumode's position and momentum variables) must in- 
clude (arbitrary) higher-order, nonlinear (non-Gaussian) 
transformations [8]. For this purpose, full universal- 
ity has been shown to be asymptotically approachable 
through infinite (but efficient) concatenation of a finite 
set of elementary unitaries, each lying in the neighbor- 
hood of the identity, and including at least one nonlinear 
gate Q- 

Secondly, when utilizing cluster states, in order to sat- 
isfy the above notion of full universality for CV, suffi- 
ciently large (potentially infinite) squeezing of the Gaus- 
sian cluster state is required, as otherwise the asymp- 
totic concatenation of elementary gate teleportations 
would accumulate an infinite amount of finite-squeezing- 
induced errors. The second issue here, the issue of finite 
squeezing, is then related with the first issue, the issue 
of full universality for CV based on infinite, elementary- 
gate concatenation. Although it has been proven that the 
squeezing per mode needed to create a universal Gaussian 
cluster state of fixed accuracy does not depend on the size 
of the cluster state (and hence on the size of the compu- 
tation it is used for) [lOj], the errors in a cluster compu- 
tation using a fixed-accuracy cluster would nonetheless 
grow arbitrarily with the length of the computation (and 
the size of the cluster) . 

In this paper, we focus on a restricted class of clus- 
ter computations, namely those realizing linear, Gaus- 
sian transformations corresponding to quadratic Hamil- 
tonians. More generally, these transformations are re- 
ferred to as linear unitary Bogoliubov (LUBO) transfer- 
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mations. In this case, it is well-known that arbitrary 
quadratic Hamiltonians can be exactly and finitely de- 
composed into elementary quantum optical elements such 
as single-mode squeezers and beam splitters [IJ E3 • A 
perfect simulation of the total Hamiltonian no longer re- 
quires an infinite concatenation of these elementary op- 
tical gates; each elementary gate no longer has to be 
weak and may even be far from the identity. These prop- 
erties greatly simplify the theoretical analysis and the 
experimental implementation of LUBO transformations 
through cluster computation over CV. As the Gaussian 
transformations play the roles of the Clifford gates for 
CV, the measurements in a Gaussian cluster computa- 
tion may all be done in parallel ('Gaussian parallelism'); 
moreover, local homodyne detections on the individual 
qumodes of the cluster are sufficient to achieve any multi- 
mode LUBO transformation Q. 

Despite these known simplifications and possibly be- 
cause of the known impossibility of full universality in the 
case of Gaussian cluster computations, so far there has 
been no explicit derivation of universal cluster states for 
Gaussian/Clifford computations which would include an 
explicit choice of homodyne measurements on a specif- 
ically shaped finite-sized cluster state realizing opera- 
tions far from the identity. It has only been shown how 
a single-mode squeezing transformation can be approx- 
imately applied to an arbitrary input state attached to 
a perfect (infinitely squeezed), linear four-mode cluster 
state [lj]. 

Here we shall give several such explicit derivations. In 
particular, we show that an arbitrary one-mode LUBO 
transformation can be perfectly achieved through an ideal 
four-mode linear cluster state. Further, we show that an 
arbitrary input state can be coupled to the cluster state 
using standard quantum teleportation [3, EH ■ Finally, 
we present a simple idea that enables one to implement 
an arbitrary multi-mode Gaussian transformation. Even 
though we will not give a provably optimal, multi-mode 
solution with regard to the size of the cluster, in our 
proposed scheme, the dependence of the cluster size is 
quadratic on the number of the input modes and this 
order coincides with the minimum order of elements re- 
quired for general multi-mode Gaussian transformations. 

As a consequence of our results, the efficient experi- 
mental implementation of any multi-mode LUBO trans- 
formation on any optical multi-mode quantum state (es- 
pecially including non-Gaussian input states) becomes 
possible using the existing optical schemes for efficient, 
deterministic creation of Gaussian cluster states fl6L fl7l . 
EH EU US- I n other words, the entire regime of multi- 
mode linear optical transformations becomes, in princi- 
ple, accessible through one fixed, offline squeezed, finite- 
sized cluster state and homodyne detections on it. 

The plan of the paper is as follows. First, in Sec. [Til 
we will give a brief introduction into cluster computa- 
tion over CV including the elementary teleportation cir- 
cuits for gate teleportation. In Sec. IIII1 we explicitly 
derive the linear four-mode cluster state and the homo- 



dyne measurement steps which allow for a realization of 
arbitrary one-mode LUBO transformations. In order to 
attach arbitrary quantum states to the cluster in an ef- 
ficient way, we show in Sec. IIVI how one may employ 
standard quantum teleportation for this purpose. An ex- 
plicit scheme for a one-mode LUBO transformation using 
teleportation-based input-cluster coupling is discussed in 
Sec. [V] Finally, before concluding in Sec. IVII1 we exam- 
ine the most general case of universal multi-mode LUBO 
transformations in Sec. IVI1 



II. ELEMENTARY GATE TELEPORTATIONS 

Before going into detail, we shall briefly review the ba- 
sic concepts of continuous-variable (CV) cluster compu- 
tation in quantum optics. We use the convention h = 1/2 
such that [x,p] — i/2 for a = x + ip and [a, a'] = 1, where 
the real and imaginary parts of an optical qumode's an- 
nihilation operator are as usual expressed by the position 
and momentum operators x and p, respectively. 

The building block of a one-mode cluster computation 
is shown in Fig. [TJ It can be considered as a generalized 
(kind of 'half') teleportation 0,|j|. First, the input state 
and an ancilla squeezed vacuum state \p = 0) are 
coupled through a CV quantum nondemolition (QND) 
interaction. A QND coupling between modes j and k is 
described by the gate exp(2 %X j X ), which is depicted in 
Fig. [TJ as a line that connects the two horizontal wires 
for each qumode. Next, the input mode is subject to 
a local measurement with a measurement basis {O 
(that is, the measured observable is p' = O^pO), where 
O is a function of only x, i.e., O = exp[if(x)]. After the 
feedforward operation Xj(s) = exp(—2ispj), which is a 
position displacement in phase space by the value of the 
measurement outcome s, the resulting output state cor- 
responds to \ip') — FO\ip), where F = cxp[— i(ir/2)a^a] 
is the Fourier transform operator. In the realistic case, 
\p = 0) will be approximated by a single-mode finitely 
squeezed state. As a result, some unwanted excess noise 
is introduced at each teleportation step of the computa- 
tion depending on the initial squeezing level. 

Arbitrary one-mode transformations can then be per- 
formed by concatenating sufficiently many elementary 
teleportation steps. Similarly, when several modes prop- 
agate through a two-dimensional cluster state (such as a 
2D lattice), QND gates can be applied to any two modes 
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FIG. 1: An elementary one-mode one-way QC gate; \ipi n ) is 
the input state; \p — 0} is a momentum eigenstate with eigen- 
value zero; p' is the measurement variable and X a correction 
displacement operator. 
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FIG. 2: (a) One-step elementary one-way QC gate (enclosed 
by a dashed line) and its cascade (the whole), (b) An equiv- 
alent circuit to (a). The circuit enclosed by a dashed line 
shows a four-mode linear cluster state, (c) and (d) are graph 
representations of (a) and (b), respectively. Each ball shows 
a single mode, while each thick line shows a QND connection. 
A dashed box in (d) is a four-mode linear cluster state which 
corresponds to that in (b). 

during the cluster computation such that universal multi- 
mode transformations become possible Jf|. 

Fig. [5] (a) shows an example of a cascade of teleporta- 
tion steps for one- mode manipulations. Every single step 
i will apply the operation FOi. Hence the general output 
state of an n-time cascaded one-mode circuit corresponds 
to 

FO n (x) . ..Fd 3 (x)Fd 2 (x)Fd 1 {x)\iP) in 

= F...Pd n {-) . . . d 3 {-x)d 2 {p)dx{x)\^) in , (i) 

As one can see, elementary unitary operations, either di- 
agonal in x or in p, are alternately performed on the input 
state. 

One important thing here is that the QND cou- 
pling exp(2ixjXk) is an element of the Clifford group 
C2, which is a group that consists of the normalizers 
of the Heisenberg-Weyl (HW) group Cx, i.e., C2 = 
{U\UCx& = Cx}. The HW group Cx is the group of 
phase space displacements, an element of which is gen- 
erally written in the form exp[2i YljiVj^j ~ CjPj) + *0] 
where rjj and £j are arbitrary real values that represent 
the size of the displacements in phase space for mode 
j, and is a global phase. The Clifford group C2 is a 
group whose generators are polynomials up to quadratic 
order in position Xj and momentum pj, i.e., the elements 
take on the general form exp[i J2j k ( a j,kXjXk+/3j,kXjPk + 

Ij.kPjPk) + i Ylji&j&i + e jPj) + #]) where otj,k, Pj,k, Jj,k, 
5j, and Cj are arbitrary real values. 

As a consequence of the discussion of the preceding 
paragraph, all the QND couplings can be applied prior to 



an actual quantum computation, while the feedforward 
operations remain simple displacements in phase space 
[Fig. Efb)]. The resulting multi-mode entangled state 
[see the dashed box in Fig.[2{b)], in which several single- 
mode squeezed states are coupled through pairwise QND 
interactions, is the resource cluster state. 

In the following, a cluster state built from 'blank' 
squeezed vacuum modes (i.e., without an input quantum 
state attached to it) shall be referred to as an "ancilla 
cluster state" . Once such a resource state has been pre- 
pared, the individual displacements of every teleporta- 
tion step can then all be postponed until the end of the 
cluster computation, as illustrated in Fig. [2fb) . How- 
ever, it does make a difference whether the desired op- 
eration Oj G C2 or Oj ^ C%. In the latter case, when 
Oj C2 for some j (corresponding to cubic or higher- 
order gates), the measurement bases of the succeeding 
((j+l)th, (j+2)th, . . . ) teleportation steps would depend 
on the outcome of measurement j. More conveniently, 
when Oj € C2 for all j , none of the chosen measurement 
bases depend on any measurement outcomes such that 
all the measurements can be performed in any order. 

Cluster states are often represented using graphs [2l| . 
as, for example, the four-mode linear chain in the dashed 
box of Fig. [2fd) where each node denotes an ancilla 
single-mode squeezed state and each link represents a 
QND coupling. Using such graphs, we can easily distin- 
guish different types of entangled cluster states. A per- 
fect cluster state can be approached in the limit of infinite 
ancilla squeezing with the resulting quantum correlations 
for all j IH, 

Pj- E f *^°> ( 2 ) 

keN(j) 

where N(j) denotes the set of all nearest neighbors to 
the j-th mode. In the limit of infinite squeezing, these 
quantum correlations among the qumodes' quadratures 
uniquely determine the corresponding graph state. The 
correlations are analogous to the generators of the stabi- 
lizer group for a qubit graph state The only differ- 
ence here is that for CV, it is more convenient to express 
the stabilizer conditions in terms of the Lie algebra, i.e., 
the generators of the HW Lie group, for which the stabi- 
lizers become 'milliners' (Ioj |. 

In the following, we restrict ourselves to unitary Gaus- 
sian transformations on n modes, which form a Clifford 
group C2 — Cl(n). The Clifford group is a semidirect 
product of the symplectic group Sp(2n, M) and the HW 
group Cx = HW(n), Cl(n) = Sp(2n,M) k HW(n). The 
group HW(n) is a homogeneous space under the adjoint 
action of Cl(n), and one can construct a group repre- 
sentation of Cl(n) on the vector space of its Lie algebra 
hw(n). Here, instead of using this particular representa- 
tion, we prefer to consider a representation isomorphic to 
the former one, but revealing a clearer physical meaning: 
the linear transformation of position x and momentum p 
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in the Heisenberg picture, 

(;)»<*■. -(s ")© + (?)' <3 > 

where x (x) and p (ft) denote the vectors of position 
and momentum operators x = (xi, ■ ■ ■ , x n ) T and p = 
(pi,...,p n ) T at the input (output), respectively. The 
2n x 2n matrix M G ( n ) = ( ^ ^ ) is a faithful representation 
of the symplectic group Sp(2n,R) with 2n 2 + n degrees 
of freedom. Here, the matrix MQ/ n ) is divided into four 
n x n matrices A, B, C, and D. The column vectors 
e, / e R™ represent displacements in phase space. The 
isotropy subgroup of this representation is a global phase 
exp(i</>) which we can ignore. The displacements will 
be omitted as well, as they can be trivially applied at 
any time during a cluster computation 0, [l3| . Note that 
Eq. ^ corresponds to an n-mode LUBO transformation, 
usually expressed in terms of annihilation and creation 
operators, a k = J2i Mlh + Bkia\ +7fe, with the j k being 
n complex parameters and the n x n matrices A and B 
chosen such that the bosonic commutators are preserved. 

III. UNIVERSAL ONE-MODE LUBO 

Let us now start with the explicit realization of an arbi- 
trary one-mode Gaussian transformation M^rn ) = ( ° b d ) , 
where ad — be = 1. In cluster computation, the el- 
ementary gate for one-mode LUBO/Gaussian transfor- 
mations is the quadratic phase gate Og(x) = exp(inx 2 ) 
[22| . where k takes on arbitrary real values, together with 
the Fourier transform. Therefore, our strategy will be to 
search for decompositions of a given LUBO transforma- 
tion into quadratic x and p phase gates. In case of the x 
phase gate, the corresponding observable to be measured 
is e~ lKX pe iKX = p + kx — g(xsin9 + pcosB), where 
g = \/l + n 2 and 9 — arctanK. In an optical implemen- 
tation, any such linear combination of x and p can be 
measured by means of homodyne detection with a suit- 
able choice of the local oscillator phase depending on the 
angle 0. 

The 2x2 matrix representation of Og(x) is O(n) = 
( 1 ? ) and that of the Fourier transform F is F = 
(? "o 1 ) = W2) where R{8) = (-f is a phase 

space rotation. Thus, the total transformation of a 
one-step one-mode teleportation gate becomes M(n) = 
FO(«) = (-«- 1 ). 

Now first we prove the following lemma. This lemma 
will also be useful later on, so we describe it in a some- 
what general form. 

Lemma: Let us combine those one-mode Clifford group 
(C2) operations that are performed before the final 
elementary gate M(n n ) into M„_i(ki, . . . , K n -i) = 

(c^Iidn-O G S P( 2 ' K ) where («i, . . . , K„_i) are the 
free parameters in the choice of the measurement bases. 
Then, together with the final step M(n n ) = FO(n n ), an 



arbitrary one-mode C2 operation is accomplished if and 
only if (iff) (a„_i,6„_i) 6 R 2 covers the whole range of 
M 2 \ {0,0}. This means that a certain property of the 
whole circuit without the last step, M„_i, determines 
whether the circuit as a whole is universal or not. 
Proof: The matrix representation after the final step can 
be written as, 

M n («i, . . . , k„) = M(k„)M„_i(«i, . . . , k„_i) 

Cfi—i 1 d n — 1 K n b n — 1\ 

s (2 i) • < 4 » 

Necessity: if (a„_i, b n -\) does not cover R 2 \ {0, 0}, then 
(c„,d„) cannot take on arbitrary values in R 2 \ {0,0}, 
thus M„(ki, . . . , k„) is not universal in Sp(2,M). 
Sufficiency: in the case of c„ = a„_i 7^ 0, a n can take on 
an arbitrary real value that is determined by c„. Now a n , 
c n = fln-i 7^ 0, and d n = 6„-i take on arbitrary values, 
and b n is automatically determined from the condition 
a-ndn — b n c n = 1, as c n ^ 0. In the case when c„ — 
a„-i = 0, we have d n = 6„_i ^ 0, and M„(ki, ...,K n ) 

has the form ( ^ n ^; b n = —d n -i—K n b n -i takes on an 

arbitrary value determined by k„, as ^ 0. Q.E.D. 

Using the above lemma, we can show that the min- 
imum number of elementary steps that is required for 
universal one-mode Gaussian transformations is four. 
Because there are three degrees of freedom (DOF) for 
Sp(2,]R), one might expect that three steps are sufficient. 
However, some measure-zero set of operations in Sp(2, M.) 
cannot be achieved with only three steps. This is ex- 
pressed by the following theorem. 

Theorem: In order to realize an arbitrary one-mode 
LUBO transformation through one-way computation 
over CV, four elementary teleportation steps, involving 
quadratic phase gates and Fourier transforms, are neces- 
sary and sufficient. 

Proof: The matrix representation for two steps is 
M(k 2 )M( Ki ) = i^X'-i) = (lit), thus when 
62 = 0, the parameter 02 cannot take on a value other 
than —1. As a consequence, M(kz)M(k2)M(ki) = 

(- K ~ + -f +K1 -T 2 2+1 ) = (Zt) cannot have d 3 = 

and 63 7^ 1; hence three elementary steps 
M(k 3 )M(k 2 )M(ki) are not universal for Sp(2,R). 

On the other hand, (03, 63) = (— 63K1 + K3, — K3K2 + 1) 
does cover the whole range R 2 \ {0,0}, as follows. The 
parameter 63 takes on an arbitrary real value indepen- 
dent of K\. In the case of 63 7^ 0, 0,3 can then take 
on an arbitrary real value that is determined by n\. In 
the case of 63 = 0, K3 = 1/k 2 takes on an arbitrary 
real value different from zero, and so does 0,3. As a 
result, using the above lemma, four elementary steps 
M(k4)M(k3)M(k2)M(ki) are (necessary and) sufficient 
for universal one-mode Gaussian operations. Q.E.D. 

We complete this discussion by presenting the explicit 
choice of parameters K\, . . . , k&. The total matrix for four 
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steps is, 



M(k 4 )M(k 3 )M(k 2 )M(ki) 

V — K3K2K1+K3+/C1 — + I 



(5) 



An arbitrary one-mode Gaussian operation represented 
by Mam = ("^) € Sp(2,K) can be decomposed into 
M(k 4 )M(k 3 )M(k 2 )M(k 1 ), as follows, 



H2 



1 - d 1 - a + bK\ , . 

K3 = C — OKI, K4 = , (OJ 



K3 



K3 



where Ki is a free parameter which should be typically 
chosen such that K3 ^ 0, unless the numerators of k 2 
and K4 in the above equations are zero, for which K3 
may become zero. One simple example is the identity 
operation Mg(i) = (oi); which corresponds to «i = 
• • • = K4 = 0. 

As those operations which are not achievable through 
a three-step computation are only a small subset of the 
whole set of Gaussian operations, one might just consider 
approximations infinitesimally close to them. However, 
in the realistic case, this is not a good strategy, because 
the squeezing of the ancilla cluster states will be finite. In 
this case, the finite-squeezing-induced excess noise grows 
arbitrarily big for three-step circuits that aim at suffi- 
ciently closely approximating otherwise unachievable op- 
erations. In the four-step case, however, such large excess 
noises are avoided, and furthermore, the extra degree of 
freedom can be exploited to minimize the excess noise. 



IV. INPUT COUPLING THROUGH 
TELEPORTATION 

For Clifford one-mode one-way quantum computation, 
it is straightforward to apply the results of the preced- 
ing section on universal one-mode LUBO transformations 
directly to the most general scenario where an arbitrary 
input quantum state is attached to the ancilla cluster 
state through QND coupling. In this general case, the 
input state may have been already processed and may 
correspond to the output of an earlier quantum computa- 
tion. A crucial question then is how to achieve this input 
coupling between a fragile input quantum state and the 
ancilla cluster state in an efficient and practical way. In 
this section, we address this issue. 

Note that there is an essential difference between the 
QND couplings for the initial ancilla squeezed states and 
those that couple an input state to the ancilla state. 
Arbitrary ancilla cluster (or graph) states can be built 
through linear optics using beam splitters and offline 
single-mode squeezed states, as has been shown theoret- 
ically [l6| and also demonstrated for some examples ex- 
perimentally [18l | . Hence, as opposed to the actual input- 
cluster coupling, the QND couplings for cluster genera- 
tion can be effectively replaced by beam splitters. Of 
course, there are situations when the input state may 
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FIG. 3: Typical diagram for input coupling through quantum 
teleportation. The input state is teleported into the cluster 
state by Bell measurements on the input modes and the end 
nodes of the cluster state. 



not be coupled to the cluster from the outside. In prin- 
ciple, an arbitrary multi-mode state can be prepared as 
a subset of modes from a larger cluster state, and, in 
this case, there is no need to prepare an independent in- 
put state before the cluster-state generation. One may 
just prepare any desired state within the cluster and then 
proceed with the quantum computation. 

Such a strategy, however, can be rather inefficient, es- 
pecially, in a one-way computation with Gaussian clus- 
ter states 23]. Furthermore, there might be situations 
in which the input coupling is necessary, for instance, 
when an unknown state has been transmitted through a 
quantum channel and is to be further processed through 
cluster computation. 

Provided that efficient QND couplings are available, 
we may just prepare the ancilla cluster state offline, and 
attach an input state to the cluster through QND cou- 
pling. However, alternatively, we may also employ a non- 
local measurement for this input coupling. A so-called 
Bell measurement, which is the two-mode measurement 
used in quantum teleportation [TH, is the prime exam- 
ple for such a nonlocal measurement. In the following, 
we will discuss this type of coupling for arbitrary input 
states through quantum teleportation. In an optical re- 
alization, an important advantage is that the Bell mea- 
surement can be easily implemented with a beam splitter 
and two homodyne detections [2~jj . 

Figure [3] shows a typical diagram of the input coupling, 
in which a two-mode input state ia attached to the cluster 
through Bell measurements on the input modes together 
with suitable 'port' modes from the cluster state. We will 
discuss only this situation, though there are many other 
possible configurations that might complicate the prob- 
lem. Quantum teleportation with CV uses an Einstein- 
Podolsky- Rosen (EPR) type state as a resource. The key 
point for our proposal is that a two-mode cluster state is 
an EPR state, up to a local Fourier transform. Thus, the 
end nodes of a cluster state can be considered as an EPR 
pair, one half of which is connected to the rest of the 
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cluster state through QND interactions. By performing 
a Bell measurement on a single-mode input state and a 
cluster end node (the unconnected half of the EPR state), 
the input state will be teleported to the connected side 
of the EPR pair, located at the edge of the cluster state. 
In the case that the input state is an n-mode (entan- 
gled) state, n independent quantum teleportations using 
n cluster end nodes would couple the input to the cluster, 
as depicted in Fig. [3l 

We describe now the usual quantum teleportation pro- 
tocol for teleporting an unknown input state into a two- 
mode ancilla cluster state (an EPR state). The quan- 
tum correlations of the two-mode ancilla cluster state 
are xi — p% — ► and x 1 — pi — > 0. We choose the linear 
beam splitter transformation for the Bell measurement as 

( t ) = 5 o,i ( % ) = ^ ( \i ) ( t ) > where ^bscript '0' 
denotes the input mode, and the primes correspond the 
the output modes of the beam splitter. The input-output 
relations for this beam splitter are 



Po 



{Po 



Pi 

■ Xi 



p'x 



l = {xi -p )/v2, 
(pi+x )/V2. 



Measuring x' and x[ is equivalent to a Bell measurement 
and leads to the standard quantum teleportation without 
any extra manipulation of the input state. 

However, by modifying the nonlocal measurement ba- 
sis compared to the Bell basis, this teleportation docs 
not only couple an input state to the cluster, but it also 
manipulates the input state correspondingly. With the 
above beam splitter coupling and subsequent homodyne 
measurements, the possible operations are Gaussian, as 
we will see below. The phases of the homodyne detec- 
tions are expressed by 9o and 9\, i.e., the observables 
x' cos 6*o -(- p' Q sin 9q and x[ cos 9\ + p[ sin 9\ will be mea- 
sured. The resulting teleportation is associated with the 
following transformation: 



cos smp_+smy+ 
cos 8— cosB- 
sin#_— sin (?+ cos 0+ 




Mi 



teJ I 



(8) 



where 9± = 9 ±9i. The standard teleportation (identity 
transfer) corresponds to the case 6>o = 9\ = 0. In the case 
6>_ = 7r/2 + nir, n 6 Z, the teleportation is not success- 
ful, because one quadrature of the input state is perfectly 
measured and the information of the orthogonal quadra- 
ture is lost; correspondingly, the elements of the matrix 
M tc i (#+,#_) go to infinity. In the following, we assume 
cos#_ > 0. For the case of cos#_ < 0, we can redefine 
9' + = 9 + + 7r and 9'_ = 9- +ir, which results in identical 
transformations, i.e., M tc \(9+,9-) = M tc i(9' + ,9'_), and 
cos 9'_ > 0. 

This seemingly complicated transformation can be in- 
tuitively understood by considering the following two 
cases separately. On one hand, in the case that the two 



local measurement bases 9q and 9\ are rotated in the 
same direction and by the same amount, i.e., 0+^0 and 
0_ = 0, we obtain a phase space rotation, 



M tel (0+,O) = 



cos 9+ sin 9 



'+ 



(9) 



On the other hand, in the case that the two local mea- 
surement bases #o and 9\ are rotated in opposite direc- 
tions by the same amount, i.e., 9 + = and ^= 0, 
squeezing will occur along the 45° direction, 

M te l(O,0_) 




R(n/4)S(r(9-))R(-ir/4) 

/coshr(#_) sinhr(0_)\ 
I sinhr(6L) coshr(0_) J ' 



(10) 



where S(r) = ( cx P( r ) oxp °_ r ^ j describes a squeezing op- 
eration, with r > corresponding to p-squeezing and 
r < corresponding to x-squeezing. The squeezing pa- 
rameter r(6-) is determined by tanhr(#_) = sin^_. In 
the case of general 9 + and , the resulting operation is 
a combination of the above two cases: 

M tel (0+,0_) 

= M tcX {9+/2, 0)M tc i(0, 9-)M te i(9+/2, 0) 

= R(-9 + /2 + n/4)S(r(9-))R(-0+/2 - tt/4). (11) 

This is a 45°-tilted squeezing operation sandwiched by 
rotations at an angle of 9 + /2. In the next section, we 
will use this result to describe a general one-mode LUBO 
transformation with teleportation-based input coupling. 



V. ONE-MODE LUBO WITH 
TELEPORTATION-BASED COUPLING 

In the case that the relative phase at the beam split- 
ter (for teleportation) may be changed arbitrarily, the 
teleportation protocol alone is sufficient to realize arbi- 
trary one-mode Gaussian operations. We shall briefly 
explain this approach which partly violates the rules of 
one-way cluster protocols, as the state manipulation de- 
pends on the choice of nonlocal measurement bases (pro- 
jections onto which require corresponding adjustments of 
the beam splitter coupling for teleportation). 

It is known that an arbitrary matrix in Sp(2, R) can 
be decomposed as [l2j : 



Mao.) = R(<h.)S(Z)R(<h). 



(12) 



The corresponding LUBO transformation of the an- 



nihilation operator a is a' 



where 
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/i = cxp[— i{4>\ + 02)]cosh£ and v = exp[— i(<f>i — 
(j>2)] sinh £ . Now the 2x2 matrix representation of 
the generalized teleportation with an extra phase ro- 
tation beforehand is M tel (9 + ,9-)R(9 in ) = R{-9 + /2 + 
TT/4)S(r(9-))R(-6 + /2-n/4 + 6 in ). As a result, an arbi- 
trary one-mode Gaussian operation can be achieved with 
the appropriate choice of and 0i n . 

For the more interesting case when we stick to the rules 
of cluster computation (i.e., we consider only the DOF of 
the local measurement bases), the relative phase at the 
beam splitter must be fixed, and so an additional two- 
step quadratic phase gate followed by Fourier transforms 
is needed for an arbitrary one-mode Gaussian operation. 
In other words, when we replace the QND coupling be- 
tween the input state and the cluster state (Fig. [2]) by a 
beam splitter interaction (Fig. the required number 
of modes of the linear ancilla cluster state remains four. 

In order to show this, we use again the lemma proven 
above. We substitute 6_ by r — arctanh sin 6_ , omit the 
subscript + in 9 + , and rewrite M tc i(9 + . 



asM t ' el (0,r), 



Ki(*,r) 



cos 9 cosh r sin 9 cosh r + sinh r 

— sin 9 cosh r + sinh r cos 9 cosh r 

or 6t 



(13) 



Let us consider the loci of (or,6r) in the R 2 plane. 
When the squeezing parameter r is fixed, the locus of 
(ot,&t) = (cos 9 cosh r, sin 9 cosh r + sinh r) is a circle, 
the center of which is (0,sinhr), intersecting the or- 
axis in points (±1,0) regardless of r (Fig. [5]). Thus, 
the set of unreachable points of (ay, 6r) in R 2 \ (0, 0) 



{(a,0)|a ^ ±1,0}. As (a T ,6 T ) in M t ', 



is N 

does not cover the whole range of R 2 \ (0, 0), using the 



telV 



IVO - 
\p-o) 2 -i- 

\P=0)s- 

\P=0)- 




\A) - 
\p=6\- 
\p=o) 2 - 

\P=0) 3 - 

\P=0)< 



BS 



-H) 



(b) 



] Disp \ 



FIG. 4: (a) Input coupling scheme through quantum tele- 
portation (dashed box) followed by two elementary gates, al- 
lowing for arbitrary one-mode Gaussian operations, (b) An 
equivalent circuit to (a). The circuit enclosed by the dashed 
box in (b) shows a four-mode linear cluster state. Measure- 
ments and feedforwards can all be put at the end of the circuit. 













(r,\ J->/1 

\a ) f 




(b) r=0j 














(c) r<0 









FIG. 5: The loci of (a T , b T ) in Eq.JT3J|. (a), (b), and (c) show 
the cases of r > 0, r = 0, and r < 0, respectively. 



lemma, we conclude that an additional elementary step, 
M(k 3 ), following M' %e y{9, r) is not enough for arbitrary 
Gaussian one-mode operations. However, teleportation- 
bascd coupling followed by an additional elementary step, 
M(n 3 )MU0,r) = (~™- ~ d ^) = ^) does 
allow for arbitrary real values of (03,63) except (0,0); 
thus, using the lemma, yet another additional step M(k4,) 
added to M(K 3 )M t ' cl ((9, r) does the trick and achieves 
Gaussian one-mode universality. 

In order to prove the above statement, we have to 
show that (03,63) covers R 2 \ (0,0). For k 3 = 0, 
the unreachable points of (03,63) in M 2 \ (0,0) are 
those of (— ct, — dr)', the corresponding set is N' = 
{(0, 6)|6 7^ ±1,0}, using the same arguments as for 
(ot,6t). Therefore, by showing that an arbitrary point 
(0, 6) G N' C M 2 \ (0, 0) is attainable for some nonzero 
«3, the proof is complete. We show this as follows: for 
a 3 = 0, cos 9 should be nonzero, and K3 = —cr/ar = 
(sin 9 cosh r — sinh r)/ cos 9 cosh r. Then 63 is calculated 
as 63 = — 1/ cos 9 cosh r, which takes on an arbitrary real 
value other than zero. Q.E.D. 

Below we give the explicit choice of the measure- 
ment bases for the implementation of a particular Gaus- 
sian operation expressed as Mqm) = ("rf) through 
teleportation-based coupling followed by two additional 
elementary steps. The two parameters 9 + and (the 
measurement bases of the teleportation coupling) are de- 
termined only from the matrix elements c and d, so that 
c sin 9+ — d cos 9+ — cos 9- — c sin . Then the other pa- 
rameters are given by k 3 = — (dcos0_+cos0+)/(sin0 + + 
sin 9- ) and K4 = — [a + (cos 9 + / cos (9_ )] /c. A solution of 
these equations is, 

1 -d 

COt 6^1 = ; — , 

2c- (1 + d)cot6» 

k 3 = c — ( 1 + d) cot 6*o , 

1 — a + b cot 9q 
K4 = 



c — d cot 9q 



(14) 
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where 8q is a free parameter, which can be utilized to 
minimize excess noises, as described above. Note that 
the problem of zero denominators in the intermediate 
expressions of K3 and K4 is avoided in the final forms for 
a suitable choice of 9 . 



VI. UNIVERSAL MULTI-MODE LUBO 

In the remainder of this paper, as a final issue, we dis- 
cuss arbitrary multi-mode Gaussian operations (general 
multi-mode LUBO transformations). We will present an 
explicit way to implement any multi-mode Gaussian op- 
eration using a finite-sized cluster state and homodyne 
measurements on it. 

The one-way two-mode entangling gate proposed pre- 
viously corresponds to a QND interaction with unit 
gain (the same gate that is used to create the ancillary, 
unweighted cluster/graph state). In order to transfer 
this gate onto a two-mode input state, the state has to 
propagate through a two-dimensional cluster state. Even 
though, in principle, sufficient for achieving universality 
with CV (when supplemented by arbitrary single-mode 
gates), the use of a single fixed-gain two- mode interaction 
gate for multi-mode transformations is rather awkward, 
as arbitrary two-mode beam splitter interactions have 
proven very powerful for multi-mode linear optics 

Here, instead of a fixed-gain interaction, we pro- 
pose another type of interaction, referred to as a three- 
mode connection gate. Its configuration is shown in 
Fig. EIa),(b). In this scheme, one ancilla mode would 
function as a kind of controller of the interaction gain. 

In Fig. [6l mode inl and mode in2 represent the input 
modes (in an arbitrary, potentially entangled two-mode 
state), while mode 3, mode a, and mode b are ancilla 
squeezed vacuum modes. Mode 3 plays the role of a 
controller of the interaction; mode a and mode b are the 
end points for the propagation of quantum information 
from mode inl and mode in2, respectively. As before, 
links between cluster nodes represent QND couplings. 

The measured variable at mode 3 is 13 = 
O t (p 3 )i 3 6(p 3 ) = x + n 3 p where 6{p 3 ) = e^ 3 ^. The 
resulting interaction is exp[i?73(a;i + x 2 ) ]■ On the other 
hand, the measurements on mode inl and mode in2 
correspond to the quadratic phase gates exp(iKix 2 ) and 
e-xjp(in 2 x 2 ), respectively, followed by Fourier transforms. 
Note that the above three operators, exp [1773(0:1 + x 2 ) 2 ], 
exp(zKix^), and exp(iK2X 2 ), a ^ commute. 

As a result, by combining these three measurements, 
an arbitrary two-mode operation is achieved whose Lie 
algebra is quadratic with regard to the position operator 
x, i.e., 0(xi,x 2 ) = exp(ciif + cjXiX 2 + c 2 x 2 ). The sub- 
sequent Fourier transform effectively swaps the roles of x 
and p. Thus, by cascading such three-mode connection 
gates, as illustrated in Fig.[5Jc), the two modes effectively 
interact subsequently with alternating quadratures x and 
p for every single step. Hence an m-time cascaded inter- 




FIG. 6: (a) Three-mode connection gate. Mode inl and mode 
inl are the input modes, while modes a, b, and 3 are ancilla 
modes. The input state may be an entangled state. Measure- 
ment results on mode 1 and mode 3 are added (electronically) 
and used to displace mode a, while measurement results on 
mode 2 and mode 3 are added (electronically) and used to 
displace mode b. (b) A graph representation of (a), (c) shows 
a three-step three-mode connection gate in graph representa- 
tion. A phase-free beam splitter can be implemented using 
this configuration. Measurements and feedforwards are omit- 
ted in (c). 



action may be written as, 

F 2 ...F 2 6 m (-)...d 2 (p)d 1 (x), 



(15) 



where F 2 is a two-mode Fourier transform, and x = 
(xi,x 2 ), p = (pi, fa). 

Note that any interaction can be suppressed by setting 
773 = such that interactions may be only applied when- 
ever they are needed for a fixed cluster state. The 4x4 
matrix representation of the connection gate is, 




1 

10 

-i?3 K 2 -nz 1 




(16) 



where O2 is a 2 x 2 zero matrix and I2 is a 2 x 2 identity 
matrix; F 2 = ( -y ~q* 2 ) is the matrix representation of the 
two-mode Fourier transform. 

To complete the discussion on arbitrary Gaussian 
multi-mode transformations, we shall use the well-known 
decomposition of multi-mode Gaussian operations, usu- 
ally referred to as Bloch-Messiah reduction [T3|. An ar- 
bitrary 77,-mode Gaussian operation G, whose DOF are 
2n 2 + n, is decomposed into the form USV, where U 
and V correspond to passive linear-optics circuits with n 2 
DOF coming from beam splitters (with some fixed phase) 
and single-mode phase shifters; S represents single-mode 
squeezers applied to each mode. 

The phase shifters and squeezers are one-mode oper- 
ations which are realizable using at most four ancilla 
modes, as discussed in detail before. Thus, provided 
an explicit implementation of a phase-free beam splitter 



with arbitrary reflectivity R is given, we can conclude 
that any multi-mode Gaussian operation is achievable 
with our specifically shaped, finite-sized cluster (where 
our implementation may be suboptimal). A decompo- 
sition of the linear-optics circuits U and V into beam 
splitters and phase shifters requires at most n(n — l)/2 
phase- free beam splitters and n(n + l)/2 phase shifters 
[111 ]. Thus, the number of ancilla modes required for 
this implementation is quadratic in the number of input 
modes n. It is now worth noting that the number of DOF 
of Sp(2n, R) is 2n 2 + n, corresponding to a minimum size 
of a cluster state for universal multi-mode Gaussian op- 
erations also quadratic with regard to n. Hence our one- 
way scheme with a total cluster state of size ~ n 2 (using 
a supply of four-mode linear subclusters and the corre- 
sponding subclusters for three-mode connection gates) 
would provide an efficient realization of universal multi- 
mode LUBO transformations. 

Finally, in order to establish the link between the three- 
mode connection gates and phase-free beam splitters, let 
us define a phase-free beam splitter with intensity reflec- 
tivity R, 



M R 2 

2 Mr 



Pi 

w 



Mr = 



R VT^R 

V1~r -Vr 



(17) 



Note that M| 



Mr 2 
2 Mr 



1 2 . We have the following relation, 



2 
1 2 



-1 2 
2 



i 2 o 2 

Mr 1 2 



= Mf. 



(18) 



The transformation Mj is achieved using a three-mode 
connection gate, choosing the three parameters ki,k 2 , 
and 773 in the following way, 



k 1 = VR- Vl-R, 
k 2 = -Vr- Vi-R, 



(19) 



Therefore, a phase-free beam splitter with an arbitrary 
reflectivity < R < 1 can be implemented through a 
three-step three-mode connection gate. This would re- 
quire in total nine ancilla modes. 



sufficient resource for universal one-mode Gaussian oper- 
ations. We have also presented a strategy for multi-mode 
Gaussian operations, where beam splitter interactions are 
used as the sole multi-mode operation. Arbitrary (phase- 
free) beam splitters can be realized in a measurement- 
based one-way scheme through so-called three-mode con- 
nection gates consuming one ancilla three-mode cluster 
per gate. Every beam splitter requires three such three- 
mode connection gates, so nine ancilla modes in total. 

Most importantly, our scheme scales quadratic with 
the number of input modes such that an ancilla clus- 
ter state of size at most quadratic in the number of 
input modes is sufficient. This scaling coincides with 
the scaling of the number of elementary optical gates 
(phase shifters, beam splitters, and squeezers) needed for 
a circuit implementation of general LUBO transforma- 
tions. We leave a possible optimization of our multi-mode 
cluster-based scheme for future research. 

Towards actual experimental demonstrations of the re- 
sults derived here, we discussed some simplifications for 
coupling arbitrary input states to a given cluster state. 
Our simplified scheme would be based on standard quan- 
tum teleportation instead of the more expensive QND 
coupling. Remarkably, eventually, the coupling QND 
gate may just be replaced by a fixed beam splitter, as 
already through our generalized teleportation scheme, it 
is possible to manipulate and process the input state to 
some extent. 

One big strength of our scheme is as follows. As it 
is well-known how to generate arbitrary cluster/graph 
states using linear optics, by employing the present 
scheme, one may now perform a general multi-mode 
LUBO transformation on an arbitrary multi-mode input 
state (including fragile non-Gaussian states) in an effi- 
cient, solely measurement-based fashion. All potentially 
inefficient, optical interactions (such as online squeezing) 
would be done beforehand offline for the resource clus- 
ter state. Although efficient multi-mode LUBO transfor- 
mations are now, in principle, accessible even for non- 
Gaussian input states, in a realistic scheme, only an ap- 
proximate, finitely squeezed ancilla cluster state could 
be used. Therefore, the resulting LUBO transformations 
would become imperfect, depending on the initial squeez- 
ing level. Apart from utilizing new experimental schemes 
with further increasing squeezing levels, one possibility 
to address the finite-squeezing issue may be in form of 
some kind of error correction such as postselection [f| or 
redundant encoding [iff]. 



VII. CONCLUSION 

In conclusion, we have described an explicit implemen- 
tation for arbitrary one-mode and multi-mode linear uni- 
tary Bogoliubov (LUBO) transformations (Gaussian op- 
erations) in the framework of one-way computation over 
continuous variables using Gaussian cluster states and 
homodyne measurements. We have shown that an an- 
cillary, linear four-mode cluster state is a necessary and 
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